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A two-component four-fermion system with equal masses, interspecies s-wave scattering length as 
and vanishing intraspecies interactions under external spherically symmetric harmonic confinement 
is considered. Using a correlated Gaussian basis set expansion approach, we determine the energies 
and various structural properties of the energetically lowest-lying gas-like state throughout the 
crossover for various ranges of the underlying two-body potential. Extrapolating to the zero-range 
limit, our numerical results show explicitly that the total energy, the trap energy as well as certain 
aspects of the pair distribution function and of the momentum distribution are related through the 
so-called integrated contact intensity I (as). Furthermore, it is shown explicitly that the total energy 
and the trap energy are related through a generalized virial theorem that accounts for a non-zero 
range. 

PACS numbers: 



I. INTRODUCTION 

Dilute equal-mass two-component Fermi gases at uni- 
tarity, i.e., gases with infinitely large interspecies s-wave 
scattering length Os and vanishing intraspecies scatter- 
ing length, show unique universal properties P, [3, H, 0, 
H,M 0, B @- Here, the term dilute means that nrg is 
much smaller than 1, where tq denotes the range of the 
underlying two-body potential and n the density of the 
system (for trapped gases, n stands for the peak den- 
sity). Intuitively, the existence of universal properties 
of the unitary two-component Fermi gas can be under- 
stood by realizing that the infinitely large s-wave scat- 
tering length as does not establish a meaningful length 
scale for the system. Correspondingly, the only length 
scale of the infinitely strongly interacting system is set 
by the density (for homogeneous systems) or the trap- 
ping potential (for inhomogeneous systems). The fact 
that the unitary two-component Fermi gas is character- 
ized by a single length scale makes it, in certain respects, 
similar to the non-interacting system. The above argu- 
ment has been formalized within a hyperspherical frame- 
work [8|. It has been shown, e.g., that the harmoni- 
cally trapped two-component Fermi gas is characterized, 
just as the non-interacting system, by ladders of states 
whose energies are spaced hy 2nuj where to de- 

notes the angular trapping frequency. The 2Huj spacing 
has been verified semi-analytically for the trapped two- 

and numerically for 
. In addition, a num- 



and three-fermion systems [111 , 
systems with up to six fermions 
ber of other universal relations have been established at 
unitarity [J, 1 i i, i, i, 0, 1 i • 

More recently, the notion of universal relations has 
been extended to dilute equal-mass two-component 
Fermi gases with finite interspecies scattering length as 
for which the two-body range vq is small compared to 
as [M [11,0, ,17, 18, IB]. While it has been established 
early on that the behaviors of these systems through the 
crossover is governed by the s-wave scattering length as 



alone (see Ref. 12(|| for a review) and that the system 
is, in this sense, universal, only recently has there been 
significant progress in relating seemingly disconnected 
quantities such as the total energy, the expectation value 
of the trapping potential, the pair distribution function 
and the momentum distribution through universal rela- 
tions that involve a quan tity termed the integrated con- 
tact intensity I{as) \l4 . [Tsl . [l6j . Physically, the inte- 
grated contact intensity measures the "number of pairs 
per unit length" . It vanishes for vanishing as , indicat- 
ing the absence of pairs, and increases as the system en- 
ters the strongly-interacting and eventually the weakly 
repulsively-interacting BEC regime, where the pairs turn 
into actual molecules [2l|, . 

This paper explicitly tests, to the best of our knowledge 
for the first time, the equivalence of four definitions of 
the integrated contact intensity /(a^) [3, [H, [l^. In do- 
ing so, we explicitly determine the integrated contact in- 
tensity for the energetically lowest-lying gas-like state of 
the four-fermion system. The spin-balanced four-fermion 
system possesses non-trivial correlations and can, at least 
in principle, be realized by loading a degenerate Fermi 
gas consisting of dimers into an optical lattice (see, e.g., 
Ref. [ill)- Once a situation with one molecule per site is 
realized, neighboring sites can be merged and molecules 
dissociated by changing the s-wave scattering length. 
This scheme would allow for the realization of an array 
of four-fermion systems with adjustable s-wave scattering 
length as- The integrated contact intensity determined 
in Sec, mil combined with the imaginary part of the two- 
body s-wave scattering 18], readily leads to a prediction 
of the atom losses due to inelastic two-body collisions, 
an experimentally observable quantity. Our results thus 
provide not only an illustration of the equivalence of var- 
ious definitions of I (as) but also pave the way for making 
quantitative predictions. 

Section |TT] reviews the definition of the integrated con- 
tact intensity I (as) for a trapped two-component Fermi 
gas interacting through a zero-range potential with s- 
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wave scattering length a^. Section IIIII illustrates the 
equivalence of various definitions of the integrated con- 
tact intensity numerically. To this end, we solve the 
Schrodinger equation for the four-fcrmion system as a 
function of the range of the two-body potential and ex- 
trapolate our solutions to the zero-range limit. Sec- 
tion IIVI applies the generalized virial theorem derived 
by Werner [l7j to our finite-range four-body energies. 
Lastly, Sec. fVl concludes. 



II. DEFINITION OF THE INTEGRATED 
CONTACT INTENSITY 

This section reviews the definitions of the integrated 
contact intensity I [as) introduced by Tan [ij, [HI, HBl- 
For concreteness, we consider a two-component equal- 
mass Fermi gas consisting of N atoms under external 
spherically symmetric harmonic confinement in which 
unlike fermions interact through a short-range two-body 
potential Vtb- For this system, the model Hamiltonian 
H can be written as 



N r 



JTi o 1 o o 



J = l k=l 



where m denotes the atom mass, lu the angular trap- 
ping frequency, iVj the number of spin- up atoms, the 
number of spin-down atoms, and fj the position vector 
of the jth atom measured with respect to the center of 
the trap. The spherically symmetric two-body poten- 
tial Vth depends on the interparticle distance rjk-, where 
"i^jk = — '^fel- Throughout, we assume that Itb is 
characterized by two length scales, the s-wave scattering 
length Os and the range tq; this assumption implies, e.g., 
that the p-wave scattering length and other higher partial 
wave scattering lengths are zero. The lengths as and ro, 
together with the harmonic oscillator length aho, where 
aho = constitute the relevant length scales of 

the Hamiltonian H. Section IIIII discusses our solutions 
to the stationary Schrodinger equation Hip = Eip, where 
the many-body wave function depends on the coordi- 
nates of all particles, i.e., ip = tpifi) ' ' ' j r^f)- Through- 
out, we indicate the dependence of the eigenenergy E on 
the s-wave scattering length and the range ro of the 
underlying two-body potential explicitly, i.e., we write 
E = E{as,rQ). 

The integrated contact intensity I {as) introduced by 
Tan [3, m, HBl applies to the Hamiltonian given in 
Eq. ID) in the limit that the two-body potential Vth is 
characterized by a vanishing range tq. Although the zero- 
range limit is, strictly speaking, not realized in nature, 
cold atom systems come arbitrarily close: The range of 
the two-body potential between two alkali atoms is set 
by the van der Waals tail, which — to leading order — falls 
off as Ce / r^jf. ; typical values of the van der Waals length 
for alkali atoms are of the order of lOOoo [l^]. The zero- 
range limit is thus approximately realized if the absolute 



value of the s-wave scattering length as and the har- 
monic oscillator length a^o are much larger than the van 
der Waals length. This scenario can be realized in cold 
atom experiments by choosing a sufficiently small trap- 
ping frequency uj and by tuning the s-wave scattering 
length through the application of an external magnetic 
field in the vicinity of a Fano-Feshbach resonance. In 
the zero-range limit, the two-body interactions impose 
the Bethe-Peierls boundary condition whenever two un- 
like particles approach each other (see, e.g., Ref. [20j). 
making even the strongly-correlated many-body problem 
amenable to analytical treatments. 

The integrated contact intensity I {as) can be defined 
through four distinct relationships [see (i) through (iv) 
below] [3 [3j HBl; as such, the integrated contact in- 
tensity establishes non-trivial connections between prop- 
erties such as the energy, the pair distribution function 
and the momentum distribution of the two-component 
Fermi gas with zero-range interactions. For concreteness, 
we restrict ourselves to the zero temperature situation, 
where expectation values are calculated for a given many- 
body eigenstate tp{fi,- ■ ■ jf^). We note, however, that 
the Tan relations can be generalized to finite tempera- 
ture [3]. 

(i) The adiabatic energy relation [3 . [isl . [3] defines the 
integrated contact intensity /(a^) through 



-^adia (^s) 



Anm dE{as,0) 
d{-a7')' 



(2) 



(ii) The virial relation [3, 13, 13] defines the integrated 
contact intensity I {as) through 

/viriai(a.) = ^^P^ [2Vt,{as,0) - E{as,0)] , (3) 

where VtT{as,ro) denotes the expectation value of the 
trapping potential, Vti = (^■mu)'^r'^j/2). We refer to 
Vtr as trap energy. 

(iii) The pair relation (3 . [isl . [3 ] defines the integrated 
contact intensity in terms of the number N^^^^ of up- 
down pairs with distances r smaller than s. 



/pair(a,s) = hm 



47r N:< 



pair 



s^Q S 



(4) 



If ^pair('") denotes the pair distribution function for the 
up-down distance r, normalized such that 



47r / Pp,i,{r)r^dr = N^N^, 



(5) 



then Eq. Q can alternatively be expressed as 

/pair (as) = lim(47r)2 Pp,i,{r) r\ (6) 

where Ppair('') is understood to be calculated for zero- 
range interactions. 
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(iv) The momentum relation [3, [T^ [TB| defines the in- 
tegrated contact intensity in terms of the number N^^^ 
of atoms with momentum k greater than K, 



Ikias) 



hm 



K- 



(7) 



This momentum relation can be reexpressed in terms of 
the momentum distribution ni{k) of the up-atoms and in 
terms of the momentum distribution ni{k) of the down 
atoms: 

(iva) If n^{k) is normahzed such that 



n^{k)(fk = N^, 



then Eq. ^ can alternatively be expressed as 

T f \ r o 2 . nao^iik) ^ 
J-k,\[0's) — hm zvr 47r ^ k , 



(8) 



(9) 

where noo.]{k) is defined through the partial wave expan- 
sion 



(10) 



and where (fc) is understood to be calculated for zero- 
range interactions. The partial wave decomposition of 
n|(/c) has been introduced to emphasize that /fej(a<;) 
only depends on the spherically symmetric component 
of n|(fc), i.e., only on k and not on k. 

(ivb) Alternatively, Eq. ^ can be rewritten in terms 
of the momentum distribution ^^(fc) of the down-atoms. 
To this end, the subscript "f in (iva) needs to replaced 
by"i". 

For N^i — Ni, we have n^{k) 



ni{k) = n{k) 
noo(fc), and consequently 
For spin-imbalanced two- 



and 71oo,t(A;) = "-ooa(^) 
h,^{as) = Ikdio-s) = Ikias)- 
component Fermi gases, in contrast, relations (iva) and 
(ivb) establish a non-trivial relationship between the mo- 
mentum distributions of the up- and down-atoms. 

The quantities /adia, -^viriai, lBmri_Ik,'[ and Ik,i have 
been predicted to be identical [Tj, |T^. To the best of 
our knowledge this has, however, not yet been illustrated 
explicitly. In fact, the determination of the integrated 
contact intensity has so far only been pursued for a few 
selected systems and limiting cases [ij, [H, [H, [H, [l^ . 
Section Hill presents explicit calculations for the energeti- 
cally lowest-lying gas-like state of the four-body fermion 
system throughout the crossover and numerically verifies 
the equivalence of the definitions (i) through (iv). 



III. INTEGRATED CONTACT INTENSITY FOR 
FOUR-FERMION SYSTEM 

This section determines the integrated contact inten- 
sity I{as) for the energetically lowest-lying gas- like state 



of the trapped four-fermion system with Nt 



as a function of the s-wave scattering length ac- 
cording to definitions (i) through (iv). Unlike for the 
three- fermion problem (see, e.g., Ref. [25]), no analyt- 
ical solutions are known for the four-fermion problem 
with zero-range interactions. Here, we solve the sta- 
tionary Schrodinger equation for the four-fermion sys- 
tem numerically u sing the correlated Gaussian (CG) ap- 
proach [H m, M, mM^ 

The CG approach[13n26|, [27|, [28l [29t expands the wave 



function in terms of Gaussian basis functions <&/ 



V'(ri,- • • ,fAr) = [-4$/(ri2,- • • ,rN~iM)] , (H) 

1=1 

where the Cfc denote expansion coefficients and 



^i{ri2, ■ ■ ■ ,rM-i,N) = exp 



N 

j<k ^^'Jfe 



(12) 



Each basis function <I>; is parameterized by N{N — l)/2 
widths dijk (one for each interparticle distance), which 
are optimized semi-stochastically following the ideas out- 
lined in Ref. [2^. In Eq. pT|) . B denotes the number 
of unsymmetrized basis functions. For the four-fermion 
system with two up- and two down-atoms, the anti- 
symmerization operator A can be conveniently written in 
terms of the permutation operator Pjk, which exchanges 
the position vectors of the jth and kih atom. If the po- 
sition vectors fi and r2 belong to the up-atoms and the 
position vectors and to the down-atoms, then A can 
be written as 



A^l-Pv 



P: 



34 



P^■,,P: 



12^34- 



(13) 



The basis functions $/ have vanishing total angular mo- 
mentum L and positive parity tt and are thus well suited 
to describe the energetically lowest-lying gas-like state 
of the four-fermion system which has L'^ = 0"*" symme- 

try [a am . 

We parametrize the two-body potential between each 
pair of up-down atoms by a simple Gaussian potential 
with range tq and depth Vb, 



Vth{r) = -Vbexp 



2rl 



(14) 



For this interaction potential, the Hamiltonian matrix el- 
ements can be calculated analytically 26]. Throughout, 
we restrict ourselves to two-body potentials that support 
either no or one s-wave two-body bound state in free- 
space. For a given range tq, the depth Vb (Vb > 0) is ad- 
justed so that the scattering length Os takes the desired 
value. While Vth leads, in general, to non-vanishing higer 
partial wave scattering lengths for finite tq, the impor- 
tance of these higher partial wave contributions decreases 
with decreasing ro and vanishes for tq = 0. To mimick 
the zero-range limit, we perform calculations for various 
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ro while keeping the s-wave scattering length fixed and 
extrapolate the quantity of interest to rg = (see below) . 
For the parameter combinations considered in this paper, 
we find that the ground state of the four-fcrmion system 
has gas- like character, i.e., that self-bound trimer and 
tetramer states are absent, in ag reement with findings 
for zero-range interactions [30|, |3l| . 

The CG approach results in an upper bound to the 
exact ground state energy. Our calculations reported be- 
low employ between B = 450 and 500 basis functions; 
we have checked that larger basis sets do not lead to a 
significant reduction of the energy. Once the expansion 
coefficients q = 1, • • • , _B) have been determined by di- 
agonalizing the generalized eigenvalue problem (the lin- 
ear dependence of the basis functions gives rise to a non- 
diagonal overlap matrix), we calculate other expectation 
values such as Vtr{as,ro), Ppair{r) and noo(fc). The ma- 
trix elements associated with these observables can be de- 
termined analytically and are readily implemented. The 
computational time for the determination of the ground 
state energy and the structural properties for a given as 
and To amounts to about an hour on a single processor 
of a state of the art desktop computer. 

Our calculations cover the weakly-attractive and 
weakly-repulsive regimes as well as the strongly- 
interacting unitary regime; in particular, we consider 41 
scattering lengths in the interval aho/fls S [—10, 10]. For 
negative scattering lengths as, we find that both the to- 
tal energy E{as,ro) and the trap energy 14r(as,''o) vary 
linearly with rg for fixed Og . For negative Os , we typically 
consider three different ranges ro, i.e., tq = 0.01, 0.03 and 
O.05aho7 and determine the zero-range quantities E(as, 0) 
and Vtr(as,0) by performing a linear two-parameter fit. 
We checked for selected scattering lengths that the in- 
clusion of additional ranges leaves the fitting parameters 
essentially unchanged. As the inverse scattering length 
aj^ is increased to positive values, the total energy starts 
depending non- linearly on the range rg. In this regime, 
we consider six different ranges rg (rg < O.OSoho) for each 
as and determine the zero-range limit E{as,Q) by per- 
forming a quadratic three-parameter fit. We find that a 
three-parameter fit provides a reliable description of the 
trap energy ]4r(as, rg) for Oho/os = to 2. For larger in- 
verse scattering lengths, Vtr(asj ''o) is significantly smaller 
than the absolute value of the total energy and has thus a 
comparatively large uncertainty. For aho/fts > 2, we de- 
termine the zero-range limit of Vtr (as, rg) by performing 
a linear two-parameter fit. While Vtr(as, rg) may depend 
non-linearly on rg in this regime, our CG data do not 
allow for a reliable determination of the functional de- 
pendence. 

As an example, circles and squares in Fig. [Ija) show 
E{as,ro) and 2Vtr(as,rg) in units of the oscillator en- 
ergy Eho, E]io — huj, obtained by the CG approach as a 
function of the range rg at unitarity. Figure UKa) shows 
that the CG quantities are well described by a linear fit 
(solid lines). Furthermore, the four-body results shown 
in Fig. [TJa) demonstrate the validity of the virial theo- 




FIG. 1: (Color online) (a) Circles and squares show E{as,ro) 
and 2Vti{as,ro), calculated using the CG approach, as a func- 
tion of the range ro of the two-body interaction potential at 
unitarity, i.e., for aho/ds ~ 0. Solid lines show linear fits to 
the CG data, (b) Circles show E{aa,ro), calculated using the 
CG approach, as a function of ro for a^^o/cLs = 5. The solid 
line shows a quadratic fit to the CG energies. For both panels, 
the basis set error of the CG data is estimated to be smaller 
than the symbol size. 

rem for zero-range interactions at unitarity 0, H, [H, , 
which states that the energy equals twice the trap energy. 
We find that the linear fit to i5(oo,rg) and 2Vi;r(oOj ^o)j 
which includes all seven data points shown in Fig. [Ija), 
predicts the same intercepts within errorbars. In partic- 
ular, we find E{oo,0) = 5.0115(5)fiiJ and 2ytr(oo,0) = 
5.0118(6)ftw, where the number in round brackets indi- 
cates the uncertainty that results from the fit alone (ex- 
cluding possible basis set errors of E and Vti themselves) . 
Our energy E(oo, 0) is slightly higher than the value of 
5.0096ftti; of Ref. [3J]; this slight difference can be at- 
tributed to a tiny basis set error. 

To illustrate the comparatively strong dependence of 
E{as,ro) on rg for positive scattering lengths. Fig. [TJb) 
shows the total energy as a function of rg for Uho/ a.s ~ 5. 
In this case, the three-parameter fit results in a zero- 
range energy E{as = aho/5,0) of — 46.93(3)/ia;, where 
the number in round brackets indicates the uncertainty 
that results from the fit alone (excluding possible ba- 
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sis set errors of the energies themselves). This uncer- 
tainty is larger than the uncertainty of Vti-{as, 0); in par- 
ticular, we find Vtr(as = aho/5,0) = 1.619{2)huj. Even 
though the basis set extrapolation error associated with 
Vtr(as, fo) might be somewhat larger than that associated 
with E{as, ro), our analysis suggests that the main uncer- 
tainty of the integrated contact intensities /adia(is) and 
-^viriai(as) m the Small positive scattering length regime 
originates in the difficulty of determining the zero-range 
limit of E(as,rQ) more reliably. 

The extrapolated zero-range quantities E{as,0) and 
Vtr(as,0) determine the integrated contact intensities 
ladmio-s) and /viriai(as) [sce definitions (i) and (ii) above]. 
To determine /adia(as)j we interpolate i?(as,0) and cal- 
culate its derivative with respect to aj^ based on this 
interpolation. Solid lines in Figs.ma)-[2I^c) show the inte- 
grated contact intensity /adia(as) as a function of the in- 
verse scattering length aj^. Figure [2] shows that /adialas) 
increases monotonically with increasing a^^. This mono- 
tonic increase reflects the increase of the two-body at- 
traction with increasing aj^. In going from the weakly- 
attractive regime (flho/as < —10) to the weakly- repulsive 
regime (flho/as > 10), the integrated contact intensity 
changes by about three orders of magnitude. For com- 
parison, dotted lines show the integrated contact inten- 
sity /viriai(os). The two contact intensities /adia(Qs) and 
^viriai(as) are nearly indistinguishable on the scales shown 
in Figs.[Ha)-[2lc), thereby lending strong support for the 
equivalence of definitions (i) and (ii). 

To quantify the agreement of /adia(as) and /viriai(as), 
Fig. [3] shows the fractional difference between these two 
quantities. Figure [3] shows that the fractional differ- 
ence is very small (less than 0.004) for negative scat- 
tering lengths where our extrapolated zero-range quanti- 
ties E(as,0) and Vti-{as,0) have the smallest uncertain- 
ties. For positive scattering lengths, the fractional differ- 
ence shows a fairly systematic deviation from zero that 
might be attributed to small systematic errors of E{as,0) 
and Vtr(as,0). The better than 1% agreement between 
^adia(as) and /viriai(as) ovcr the entire scattering length 
range considered lends strong numerical support for the 
prediction that Eqs. ^ and ^ constitute equivalent def- 
initions of the integrated contact intensity /(os). 

Next, we discuss the determination of the integrated 
contact intensity from the pair distribution functions 
^pair('') for the up-down distance. Figure 2] shows exam- 
plary scaled pair distribution functions for tq — O.OOSaho 
and various s-wave scattering lengths (see also 
Ref. [29]). In particular. Fig. gja) shows (47r)^Ppair(r)r^ 
for negative Ug while Figs.lD^b) andlD^c) show respectively 
the small r and the large r regions of (47r)^Ppair(f )?'^ for 
positive tts ■ The overall behavior of the pair distribution 
functions can be described as follows: (i) The amplitude 
at large r decreases with increasing aj^ , reflecting the 
fact that the four-fermion system becomes more compact 
with increasing attraction, (ii) As aj^ increases from 
— lOa^^ to (and to even larger values), the scaled pair 
distribution functions develop a two peak structure that 
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FIG. 2: (Color online) Integrated contact intensities 7adia(as), 
7viriai(as), -^'pair(as) and Ik{a,s) as a function of the inverse s- 
wave scattering length aj^. Solid and dotted lines show the 
integrated contact intensities /adia(«s) and 7viriai(as) calcu- 
lated according to Eqs. ((2]) and Q, respectively, using the 
extrapolated zero-range quantities as input (the two data sets 
are nearly indistinguishable on the scale shown and compared 
in more detail in Fig. [Sjl . Circles and squares show the inte- 
grated contact intensities /pair(as) and Ik{as), respectively, 
determined from the pair distribution functions and the mo- 
mentum distributions for ro = 0.005aho (see text for details). 
Panel (a) shows the entire crossover region on a logscale, while 
panels (b) and (c) show the negative and positive scattering 
length regions on a linear scale. 



reflects the formation of pairs. The peak at smaller inter- 
particle up-down distances r corresponds to the forma- 
tion of a pair while the peak at larger r values (r « Icho 
to 1.5aho) reflects the fact that the second up-atom and 
the second down-atom are pushed away from the first pair 
due to the Pauli exclusion principle (i.e., the formation 
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FIG. 3: (Color online) Fractional difference [/viriai(as) — 
Ia.dis.{cLs)]/ Iiidia.{a.s) as a function of the inverse s-wave scat- 
tering length aj^ . 



of self-bound trimers and tetramers is prohibited), (iii) 
The scaled pair distribution functions vanish at r = for 
all s-wave scattering lengths; for numerical reasons, the 
first point of Ppaii ('') is not calculated at r = but at 
r = 0.005aho- The vanishing of Ppah {r)r^ at r = is 
accompanied by a sharp drop of Ppaii{r)r^ at r values of 
the order of a few times the range ro- The scaled pair 
distribution functions behave universally when r is larger 
than a few times the range rp; for smaller r, the pair dis- 
tribution functions acquire non-universal behavior. 

To illustrate the universal, range independent behav- 
ior of the scaled pair distribution functions. Fig. [5] shows 
the quantity (47r)^Ppair('')?'^ for a number of different ro 
but fixed s-wave scattering length, i.e., for aho/is = —10 
(corresponding to — — O.laho)- For this scattering 
length, the condition vq <^ ja^l is approximately full- 
filled if To < 0.02aho. In agreement with this condi- 
tion, the universal part of the pair distribution functions 
Ppair(''), i.e., the part where Ppair(?') is independent of 
ro, extends to smaller r values with decreasing tq. Fig- 
ure [5] shows that the scaled pair distribution functions 
for small ro, ro < 0.02aho, behave approximately linearly 
in the regime 3ro ^ r < 0.05 — O.laho; in this regime, 
essentially no dependence on ro is visible. The behavior 
of Pp&iri^)'!"^ , illustrated exemplarily for = —O.laho 
in Fig. [5l suggests that the integrated contact intensity 
^pair(as) can be determined readily by employing a linear 
two parameter fit to the small r-regime where the uni- 
versal part of (47r)^Ppair(r)r^ varies linearly and where 
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FIG. 4; (Color online) Scaled pair distribution functions 
as a function of r for ro — O.OOSaho- Panel (a) 
shows the scaled pair distribution functions for aho/as ~ —10, 
—5, —2.5, —1.5, —1, —0.5 and (from bottom to top at 
r = 0.25aho). Panels (b) and (c) show the scaled pair distri- 
bution functions for r < O.Soho and r > O.Saho, respectively, 
for flho/tts = 0.5, 1, 1.5, 2.5, 5 and 10 [from (b) bottom to top 
at r = 0.05aho and (c) top to bottom at r = 2aho]. Note the 
different scales of the axes in panels (a) through (c) . 



that varies linearly extends over a comparatively large 
r range and where the slope is comparatively shallow. 
For small \as\, the extrapolated integrated contact in- 
tensity /pair(as) shows & notable dependence on the fit- 
ting range employed. The integrated contact intensity 
Ppair('') is calculated for a finite but sufficiently small /pair (is), determined from the pair distribution func- 



ro. Alternatively, one might consider extrapolating the 
pair distribution functions themselves to ro = 0; this ap- 
proach is not pursued in this work. 

Figures [n{a)-[n]Jc) illustrate our determination of 
^pair(as) for aho/as = —5, and 5. We find that 
Ip^^il■{o■s)^ which is determined by the r = value of 
the fit (i.e., the intercept), can be determined most reli- 
ably for large \as\ where the universal part of Ppair(^)''^ 



tions for ro = 0.005aho, is shown by circles in Fig. [2l 
While /pair(as) compares very favorably with /adia(as) 
and /viriai(as) in the strongly-interacting regime, /pair(as) 
is notably smaller than /adia(as) and /viriai(as) in the 
weakly-attractive and weakly-repulsive regimes. It is 
clear from Figs. EH a) through [S]^c) that /pair(as) would 
be somewhat larger in the weakly-interacting regimes if 
the upper fitting limit was reduced somewhat (see also 
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r/ a, 



ho 



FIG. 5: (Color online) Scaled pair distribution functions 
(47r)^Ppair(T')''^ as a function of r, r G [0,0.15aho], for 
aho/cis = —10 and various ro. At small r, the scaled 
pair distribution functions correspond from bottom to top to 
ro = 0.05, 0.03, 0.02, 0.01, 0.005, 0.0025, and O.OOlaho. Note 
that the solid and dashed lines for the two smallest ranges, 
i.e., for ro = 0.001 and 0.0025aho, are nearly indistinguishable 
on the scale shown. 



the caption of Fig. [S]). This would bring /pair(as) into 
even better agreement with /adia(as) and /viriai(as) in the 
weakly-interacting regimes while leaving the behavior in 
the strongly-interacting regime essentially unchanged. 

Figures[7I^a) and[7fb) show the lowest partial wave pro- 
jection noo(fc) of the momentum distribution for fixed ro, 
i.e., for ro = O.OOSaho, and selected negative and positive 
scattering lengths. For weak attraction, the lowest par- 
tial wave projection noo(fc) is approximately Gaussian. 
As aj^ increases, noo(fc) develops a tail at large k while 
the small k amplitude decreases. The tail at large k in- 
dicates that the system is characterized by increasingly 
small length scales, reflecting the formation of pairs. 

To determine the integrated contact intensity Ikia-s) 
from the lowest partial wave projection noo{k), we scale 
noo(fc) by k^ and plot the resulting quantity as a function 
of fc-i [see Eq. ^ and Figs. [Ha) and[8tb)]. Similarly to 
our analysis of the scaled pair distribution functions, we 
extract the integrated contact intensity /fc(as) by employ- 
ing a linear two-parameter fit to the linear small fc~^ re- 



gion of 27r^47rnoo(fc)/\/47r, where rtoo(fc) is calculated for 
a small but finite ro. Figures [5]Ja) and[5]Jb) illustrate the 
procedure for aho/os = —5 and 0, respectively. The re- 
sulting integrated contact intensities Ik{cLs), determined 
from the momentum distributions for ro = O.OOSaho, 
are shown as a function of the inverse s-wave scatter- 
ing length aj^ by squares in Fig. [21 The agreement be- 
tween /fc(as) and the integrated contact intensities de- 
termined based on definitions (i) and (ii) is very good 
for aj^ < 5. For larger inverse scattering lengths, no- 
table deviations are visible. These deviations can, as in 
the case of /pair(as), be traced back to the fitting range 




0.5 1 1.5 



r / a, 



ho 



"1 — I — r 
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20- 
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0.5 1 1.5 



300 
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3: 



0.5 1 1.5 



FIG. 6: (Color online) Solid lines show the scaled pair dis- 
tribution functions (47r)^Ppair(r)r^ as a function of r for 
ro — O.OOSttho and (a) aho/is = —5, (b) a^^o/as = and 
(c) aho/cia = 5. The dotted lines are obtained by performing 
a linear two-parameter fit to the small r-region of the scaled 
pair distribution functions; the fit includes r-values in the 
range [0.015aho, O.lOho]. 



employed to extract Ik{as)- Despite these deviations. 
Figs. [2ta)-[2ljc) convincingly illustrate the equivalence of 
definitions (i) through (iv) of the integrated contact in- 
tensity. 



IV. GENERALIZED VIRIAL THEOREM 

Our total energies E{as,ro) and trap energies 
Vtr(as, ''o) for the energetically lowest-lying gas-like state 
of the four-fermion system can be readily combined 
to verify the generalized virial theorem derived by 
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FIG. 7: (Color online) Scaled projected momentum distri 
bution noo(fe)/\/47r as a function of k for ro — 0.005aho 
Panel (a) shows the scaled projected momentum distributions 
for ai^o/as — —10, —5, —2.5, —1.5, —1, —0.5 and (from top 
to bottom at fc = 0.5aj~^^). Panel (b) shows the scaled pro- 
jected momentum distributions for aho/cis = 0.5, 1, 1.5, 2.5, 5 
and 10 (from top to bottom at fc = lUho^)- Note the different 
axes in panels (a) and (b). 

Werner [l7j . This theorem applies not only to infinitely 
large s-wave scattering lengths but also to finite scatter- 
ing lengths and accounts for finite range corrections, 

E{as,ra) = 

2Vty:{as,ro) - -as -ro . (15) 

2 oas 2 oro 

Here, it is assumed that the underlying two-body poten- 
tial Vth depends on only two lengths, the s-wave scat- 
tering length as and the range ro. The second term 
on the right hand side of Eq. (fT5|) can be rewritten as 
a'^^[dE{as,ro)/da's'^]/2, which shows that it vanishes at 
unitarity. In the zero-range limit (i.e., for ro = 0), 
Eq. thus reduces to the "usual" virial theorem 

E{as,0) = 2Vtr(as,0) 0, i, 113, Ull , which has been dis- 
cussed in the context of Fig. [IJa). 

Symbols in Fig. [5] show the fractional difference be- 
tween the four-fermion energy -Eviriai(as, ''o) calculated 
according to the right hand side of Eq. and the en- 
ergy calculated by the CG approach as a function of the 
inverse scattering length aj^. For negative scattering 
lengths, the virial theorem and CG energies agree to bet- 




FIG. 8; (Color online) The solid lines show the scaled pro- 
jected momentum distribution 27r'^ 4% tiqo (k) / \/47r as a func- 
tion of k~^ for ro — 0.005aho and (a) aho/as = —5 and (b) 
tSho/is — 0. The dotted lines are obtained by performing a 
linear two-parameter fit to the small fc~^-region of the scaled 
projected momentum distribution; the fit includes fc~^-values 
in the ranges (a) [0.05aho, O.loho] and (b) [O.OSaho, O.loho]. 

ter than 0.025% for the parameter combinations consid- 
ered, thus confirming the generalized virial theorem with 
high accuracy. For positive scattering lengths, our results 
confirm the generalized virial theorem with 2% accuracy. 
We note that the validity of the generalized virial the- 
orem has very recently also been investigated for large 
two-component Fermi gases by means of a Monte Carlo 
approach (ssj . 

V. CONCLUSIONS 

This paper investigates universal properties of di- 
lute equal-mass two-component s-wave interacting Fermi 
gases. In particular, our calculations are performed for 
the energetically lowest-lying gas-like state of the trapped 
four-fermion system without spin-imbalance. We deter- 
mined the energies and various other properties as func- 
tions of the scattering length and the range of the un- 
derlying two-body potential. Analyzing these quantities, 
we were able to explicitly demonstrate the equivalence of 
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FIG. 9: (Color online) Fractional difference [-Bviriai(fis, '"o) — 
E{as,ro)]/E{a3,ro) as a function of the inverse s-wave scat- 
tering length aj^ . Pluses, circles, squares and diamonds 
show the fractional difference for ro = 0.005,0.01,0.03 and 
0.05aho (for negative as, only values for the largest three ro 
are shown; for positive as, only values for the smallest three 
ro are shown). For negative fls, the absolute value of the 
fractional difference is smaller than 0.00025 for all ro shown. 



four distinctly different definitions of the integrated con- 
tact intensity I {as) [13, [H, Ell and of a generalized virial 
theorem [17] with high accuracy. In addition, we per- 
formed selected calculations for the energetically lowest- 
lying gas-like L'^ = 0+ state of the three-fermion system 
with (N-fjNi) = (2,1) and the five-fermion system with 
(iV|,iV|) = (3,2). These calculations confirm our find- 
ings discussed in Sec. lllll for the four-fcrmion system. Wc 



have outlined in detail how the integrated contact inten- 
sities /pair(as) ELud Ikio-s) Can be extracted from our nu- 
merical data. While we have not commented on how to 
best measure the integrated contact intensity experimen- 
tally, our outlined analysis provides some guidance as to 
the demands on the accuracy of experimental data if the 
integrated contact intensity is to be measured with high 
accuracy. At the same time, we have shown that the in- 
tegrated contact intensity changes by multiple orders of 
magnitude through the crossover, making a qualitative 
measurement of /(a^) appear quite plausable. 

The explicit verification of the universal relations away 
from unitarity paves the way for developing the con- 
cept of the integrated contact intensity further. For 
example, are there other properties of the many-body 
two-cornpqnent Fermi gas besides those already dis- 
cussed [3, [H, [ill that are related through the inte- 
grated contact intensity? Can the concept of the inte- 
grated contact intensity be extended to low-dimensional 
systems? Furthermore, can the concept of the inte- 
grated contact intensity be extended to relate various 
properties of unequal-mass two-component Fermi gases 
with large mass ratio or Bose gases? For these systems, 
a three-body parameter that is associated with Efimov 
physics [HjIsJ arises and it has been speculated JJ,] that 
this three-body parameter could somehow be related to 
an integrated contact intensity that, physically speaking, 
accounts for the formation of trimers. 
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